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ABSTRACT 

The shock model has successfully explained the observed behaviors of afterglows from 
long gamma-ray bursts (GRBs). Here we use it to investigate the so-called early af- 
terglows from short GRBs, which arises from blast waves that are not decelerated 
considerably by their surrounding medium. We consider a nearby medium loaded 
with e ± pairs (Beloborodov 2002). The temporal behaviors show first a soft-to-hard 
spectral evolution, from the optical to hard X-ray, and then a usual hard-to-soft evo- 
lution after the blast waves begin to decelerate. The light curves show variability, and 
consist of two peaks. The first peak, due to the pair effect, can be observed in the 
X-ray, though too faint and too short in the optical. The second peak will be eas- 
ily detected by Swift. We show that detections of the double-peak structure in the 
light curves of early afterglows are very helpful to determine all the shock parame- 
ters of short GRBs, including both the parameters of the relativistic source and the 
surroundings. Besides, from the requirement that the forward-shock emission in short 
GRBs should be below the BATSE detection threshold, we give a strong constraint 
on the shock model parameters. In particular, the initial Lorentz factor of the source 
is limited to be no more than ~ 10 3 , and the ambient medium density is inferred to 
be low, n lO^ 1 cm~ 3 . 

Key words: gamma-rays: bursts — radiation mechanisms: nonthermal — relativity 
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INTRODUCTION 



It is recognized that gamma-ray bursts (GRBs) may be di- 
vided into at least two classes: one third of the bursts with 
short duration (^ 2 s) and hard spectra, and the other two 
third with long duration (<; 2 s) and soft spectra (Kouve- 
liotou et al. 1993; Dezalay et al. 1996; Paciesas et al. 2003). 
The detections of afterglows from long/soft GRBs and then 
their redshift measurements have revealed their cosmologi- 
cal origin (see van Paradijs et al. 2000 for a review). Their 
afterglows are widely believed to come from a blast wave 
driven by a relativistic ejecta into an ambient medium (see 
reviews of Cheng & Lu [2001] and Meszaros [2002]). Unfor- 
tunately, it is impossible so far for observations to system- 
atically follow short GRBs at longer wavelengths. The effort 
of searching transient afterglow emission from short /hard 



* lizhuo@mail.ihep.ac.cn 
t daizigao@publicl.ptt.js.cn 
t tlu@nju.edu.cn 



GRB usually yields only some upper limits (e.g. Kehoe et 
al. 2001; Hurly et al. 2002; Gorosabel et al. 2002; Klotz, 
Boer & Atteia 2002). The difficulty for detection of short 
GRB afterglow is mainly due to the poor prompt localiza- 
tion by current satellites for these bursts. This problem is 
waiting for the upcoming Swift satellite to resolve. Lazzati, 
Remirez-Ruiz & Ghisellini (2001) report the discovery of a 
~ 30 s delayed, transient and fading hard X-ray emission in 
the BATSE light curves of a sample of short GRBs, the soft 
power-law spectrum and the time-evolution are consistent 
with predicted by the afterglow model. 

Based on the widely accepted blast wave model, 
Panaitescu et al. (2001) studied the long-term afterglows 
of short GRBs coming from the blast waves. In this pa- 
per, we focus on the investigation of early afterglow emis- 
sion, which arises from the blast wave before it transits to 
the self-similar evolution in the ambient medium (Bland- 
ford & McKee 1976). We consider pair loading in the exter- 
nal medium, which is caused by the collision between the 
outgoing gamma-rays and the scattered photons off the ex- 
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ternal medium (Madau & Thompson 2000; Thompson & 
Madau 2000; Dermer & Bottcher 2000; Madau, Blandford 
& Rees 2000; Meszaros , Ramirez-Ruiz & Rees 2001; Be- 
loborodov 2002; Ramirez- Ruiz, MacFadyen & Lazzati 2002). 
The pairs will affect the behavior of early afterglows. As the 
short GRBs have ~ 20 times less fluence than long GRBs 
(Mukherjee et al. 1998), the kinetic energy of short GRBs 
must be ~ 20 times less than long GRBs too, provided that 
the efficiencies for producing gamma-rays are the same for 
both classes (Panaitescu & Kumar 2001). We take the typi- 
cal kinetic energy of short GRBs as 10 52 ergs here. Further- 
more we assume that the shocks in short GRBs have param- 
eters similar to those of long GRBs, except for the ambient 
density which is believed to be lower if short GRBs origi- 
nate from the compact binary mergers (Eichler et al. 1989; 
Narayan, Paczyhski & Piran 1992). Later on we will show 
that low density for a short GRB is required (eq. (17) and 
discussions below). In section 2 we discuss the hydrodynam- 
ics of short GRBs, and in section 3, pair loading in the ex- 
ternal medium. An early afterglow from the blast wave is 
derived in section 4. Section 5 gives conclusions and obser- 
vational implications. 



2 HYDRODYNAMICS OF SHORT GRBS 

A GRB itself is believed to come from internal shocks which 
are due to different Lorentz factors of shells within the ejecta 
(Rees & Meszaros 1994). After producing GRB the ejecta 
cools down rapidly and may be considered as a cold shell. 
The interaction between the outgoing shell and ambient 
medium leads to two shocks: a forward shock propagating 
into the medium, a reverse shock sweeping up the ejecta 
matter, and a contact discontinuity separating the shocked 
ejecta matter and the shocked medium. So the kinetic en- 
ergy of the ejecta can be dissipated into the internal energy 
of the medium by the forward shock and into the internal 
energy of the ejecta matter by the reverse shock. According 
to Sari (1997), there are two time scales. One is relevant to 
the forward shock, at which the shell reaches an decelera- 
tion radius where the shell has given the medium an energy 
comparable to its initial energy, 
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where Ek = W 52 Ek,52 ergs and r\ = 300?73oo are the fireball 
kinetic energy and initial Lorentz factor, n = 0.01n_2 cm~ 3 
is the particle density of the ambient medium, and z is the 
source's redshift. The other is relevant to the reverse shock, 
at which the reverse shock accelerates to become relativistic. 
The ratio between the two time scales is defined as 
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where A is the shell width (in observer frame) of the ejecta. 
In the internal-shock model the shell width is A = cT = 
3 x 10 9 T_i cm, with T = 0.1T_i s the duration of GRB. 
Eq. (2) shows that £ is not sensitive to Ek and n, and only 
somewhat dependent on r\ which is not accepted to be quite 



larger than 10 3 (implied from eq.[17] below, and also im- 
plied from other aspects of GRBs, e.g., Lazzaiti, Ghisellini & 
Celotti 1999; Derishev, Kocharovsky & Kocharovsky 2001). 
Thus, for short GRBs, we usually have £ > 1. In this case, 
the reverse shock is initially Newtonian and becomes mildly 
relativistic when it crosses the shell at idee- Consequently 
the shocked medium has most of the initial energy, and the 
forward shock goes into the self-similar Blandford-McKce 
(1976) evolution. 



3 PAIR LOADING IN GRB MEDIUM 

The GRB itself from internal shocks is emitted early, pre- 
ceding the development of the blast wave. The gamma-ray 
front interacts with the ambient medium, leading to two 
processes: Compton scattering and 7 — 7 absorption of the 
scattered photons. As a result the medium is loaded with 
pairs within a loading radius i?i oa d = 5 x W 15 E^ 2 2 cm, with 
E-y — 10 5 -E 7 ,52 ergs the isotropically explosive energy in 
gamma-rays (Beloborodov 2002). Approximately 10 3 pairs 
per ambient electron can be created when conditions are 
right, but usually it is much less, fo = N±/N(R\ 0!l<1 ) = 
10 2 /o, 2 (Beloborodov 2002). Therefore, the mass of e pairs 
ahead of the blast wave is neglected, fo < m p /m e , it dose 
not affect the dynamics of the blast wave. Besides, the pairs 
may be pre-accelerated by the gamma-ray front, but the pair 
energy does not exceeds the ejecta kinetic energy. Provided 
that the medium density is low and the deceleration occurs 
outside the pre-accelerated radius i? acc (Beloborodov 2002) 
which is smaller than -Ri oa d, the deceleration time (eq.[l]) 
will be not affected (Note however that as shown by Be- 
loborodov [2002], for dense enough medium tdcc changes 
whenever deceleration occurs in a pre-accelerated medium, 
i.e., i?dcc < Ra.cc)- Typically, the deceleration time is longer 
than the one at which the blast wave approaches i?i oa d, 
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and the one at which the blast wave crosses a radius, Rf — 
/d^ 3 ^?ioad where the number ratio (/) of the pair to ambient 
electron number drops to / = 1, 
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Thus, for short GRBs we have the order, ^ oad <tf < t dcc . 

Here, when introducing radius Rf, we have assumed 
mixing of particles in the blast wave that allows the newly 
added post-shock particles to share energy with earlier in- 
jected pairs. Before the reverse shock crosses the ejecta and 
vanishes at tdcc, the existence of the contact discontinu- 
ity prevents the earlier pairs from being far downstream 
the forward shock front. The total shocked mediums are 
compressed between the contact discontinuity and the for- 
ward shock front. Furthermore, the coupling of leptons with 
baryons may take place in presence of even weak magnetic 
fields (e.g. Madau & Thompson, 2000; Meszaros, Ramirez- 
Ruiz & Rees 2001), therefore the total particles are possible 
to be mixing, allowing continuous transmission of energy 



© 2003 RAS, MNRAS 000, 000-000 



Temporal variability in early afterglows of short gamma-ray bursts 3 



from baryons to leptons. We will take the mixing hypothesis 
in the following. 



4 EARLY AFTERGLOWS OF SHORT GRBS 

We now derive the temporal property of early afterglows 
from forward shocks of short GRBs. We consider the source 
as an isotropic explosion even though it may have jet geom- 
etry, because the jet effect is not important at early times 
when the jet open angle is larger than ~ \jr\. 



4.1 Phase ti oa d <t <tf 

We begin with the blast wave having swept up all the pro- 
duced pairs at J?i oa d. The pairs will modify the usual prop- 
erty of the afterglow, since the same energy will be shared 
by much more leptons. Furthermore the pairs will increase 
the radiation efficiency significantly. With the mixing hy- 
pothesis, the comoving-frame random lepton Lorentz factor 



The synchrotron spectrum from leptons distributed as 
dN iep ,/d^/ e oc y~ p (7e > j m ) is a broken power-law with 
break frequencies v m and v c : F v oc u 1 ^ 3 at v < v v = 
mm(v m ,v c ); F v oc v~ x l 2 for v c < v < v m or F„ oc j, - ^- 1 )/ 2 
for v m < v < v c ; and F v oc v~ v ^ 2 at v > max(i/ m ,i/ c ) 
(Sari, Piran & Narayan 1998). Here we neglect the syn- 
chrotron self- absorption which is only important at longer 
wavelengths, e.g., radio or IR. 

Since / oc N^ 1 oc R~ 3 oc t~ 3 , eq.(6) implies that the 
peak frequency rapidly increases, as v m oc t 6 , from the op- 
tical to the hard X-ray band eventually (eq.[13]). Thus, we 
have the scaling laws for ti oad <t<tf, 

F p = const., Umtxt 6 , v c ctt~ 2 (tioad < t < tf). (10) 

The afterglow shows a soft-to-hard spectral evolution during 
this phase. Observed at a fixed frequency, v ob , between the 
optical and hard X-ray, the light curve will show a rapidly 
increase, F v oc £ 3 ( p_1 ) , and then a sharp decreasing, F v oc 
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after v m crosses v ob at 
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7m = n , Js teV- (5) 

(1 + f)m e 

Here, / = N±/N e , and the energy density in leptons and 
magnetic field — behind the shock are usually parameter- 
ized by the fractions e e = 0.1e e ,-i and es = 0.01es,-2 of 
the total internal energy density (ifnmpC 2 ), respectively. 
For / > 1 at tioad < t < tf, this Lorentz factor is a factor 
(1 + /) « / lower than usual case, and the corresponding 
synchrotron frequency is therefore 



i o inl4 2 1/2 4 
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which is in the optical band if / = fo = 10 2 /o,2 at t\ oa .d, as 
opposed to the hard X-ray of the usual case. Now the pair 
number dominates the ambient electron's number, the total 
lepton number is Mep. ~ |7ri?f oad n/o. The peak spectral 
power (in comoving frame) per lepton is P„, m ax = 1.4 x 
10~ 22 P crgss _1 Hz _1 . We then have the afterglow peak flux 



F p = N lep .r}P v . 
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where di = 10 di^s is the GRB's luminosity distance. To 
calculate the synchrotron spectrum, we still need to know 
the cooling frequency that is corresponding to those leptons 
which cool by synchrotron/inverse-Compton radiation in a 
dynamical time t, i.e. 

„r„ 1n 20 -3/2 -4 -3/2.-2 / 1 + Y\ 2 ( 1 + z \ u / oN 

v c = 2.5x10 e Bi i 2 »?3oo"'-2 t [-^—) [~2~ J Hz >( 8 ) 

where Y is the Compton parameter. According to 
Panaitescu& Kumar (2000), Y = \{[li.e e /e B ) + 1] 1/2 - 1} 
1. Now for observer's time t = ti oa d, 



^c(ftoad) = 3.4 x W 20 e-^ 2 2 E-l 2 
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4.2 Phase t f <t < t dcc 

outside Rf, we have / < 1, implying that pair effect is neg- 
ligible. The afterglow property then approaches the usual 
case, where 

Fp oc N e oc t [i , v m — const., v c oc t~ 2 (tf <t< t dcc ).(12) 
In details, 

v m = 1.8 x 10 1S e 2 i _ie B / i i 2 £fc,52r;3oo«- / 2 2 {~^~) Hz ( 13 ) 

is a constant and should be in the hard X-ray band. The 
cooling frequency continues to decrease (eq.[8]) to 

^c(tdec) =5.0xl0 1 ^ s 3 i 2 2 i5-jf % 4 / 3 n:p (^y 1 Hz(14) 

at t = tdcc- Note that v m > ^c(idcc), implying that v c has 
crossed v m at a certain moment t cm after which the spec- 
trum becomes peaking at v c , which is in X-rays. Due to 
ambient electrons picked up, the peak flux increases rapidly 
to 

(15) 



Fp(tde C ) = 580e 1 J 2 _ 2 E k ,52Ti_ 2 u, 
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at t = tdcc- If observing at a fixed sub-keV frequency, wc 
can see in this phase the light curve climbing up again. 

Now a constraint on short GRBs arises from the re- 
quirement that the flux in sub-MeV should not exceed the 
BATSE detection threshold. Otherwise, as t dcc > 2 s, the 
burst is not short any more. With n = 0.01 and other pa- 
rameters in their typical value, we obtain the flux given by 

$(MeV) ~ 2v m F Vm = 2{v m v c ) 1/2 F p = 1.1 x 10" 8 

ci2/3 8/3 1/3.-2 -2 -1 

Xe e,-i£ fc ; 52 77 3 ^n-- 2 <i (i28 ergs cm s ■ ( i6 ) 

We set that the BATSE threshold is 1 x 10~ 8 ergs cm~ 2 s~\ 
leading to a constraint on the "short GRB" parameters of 



p 2/3 8/3 1/3,-2 1 
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(17) 
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Note that the most stringent constraint is on r), which is 
not allowed to be too large, i.e. r\ 10 3 . A lower limit 
to rj arises from the requirement that during the prompt 
sub-MeV burst, the optical depth due to scattering off fire- 
ball electrons, tt, = orNb/^nR^, should be less than unity 
(Rees & Meszaros f994), with N b — E k /rjm p c 2 the fireball 
baryon number, R-, < r] 2 c5t the radius at which the fire- 
ball kinetic energy dissipated to gamma rays, and 8t the 
shortest time scale of rapid variability in the GRB profile. 
This leads to r\ > SSQE^f^StZ^ 5 ', thus the r\ value taken 
in eq. (17) is to the lower limit. If the other parameters are 
fixed to their typical values, the ambient density for short 
GRBs is limited to n ^ 0.01 cm -3 (eq. [17]), consistent with 
the clean-environment hypothesis to short GRB models of 
compact binary systems, e.g. Eichler et al. (1989); Narayan, 
Paczyhski & Piran (1992). 

4.3 Phase t > t dcc 

In this phase, the blast wave begins to decelerate consid- 
erably. If the electrons obtain a significant fraction of total 
energy, e e ~ 1, the blast wave will evolve in the radiative 
regime, since all the electrons are fast cooling, with v c < v m . 
The light curve is somewhat complicated in this case with 
light-curve index related to e e (Bottcher & Dermer 2000; Li, 
Dai & Lu 2002). For the typical value e e = 0.1, we can safely 
consider a adiabatic blast wave, so the well know scaling laws 
are: 

F p = const., v m oc t~ 3/2 , f c k i _1/2 (t > tdec), (18) 

where F p is given by eq.(15). The afterglow spectrum shows 
the usual hard-to-soft evolution after tdec- Observed at a cer- 
tain frequency v \, between the optical and keV band, when 
v m or v c crosses u oh , whichever the first, the observed flux 
reaches a peak with F b = F p ~ 580/iJy. It is a magnitude 
~ 15.6 if observed in the optical. Thus, there is another peak 
in the light curve other than the first one in the phase t <tf. 
Furthermore, this second peak is much stronger than the 
first one. Lazzati, Remirez-Ruiz & Ghisellini (2001) claim 
to have detected such a delayed hard X-ray peak. 

In Fig. 1 we show the light curves at two bands, the 
optical and the X-ray, also labelled in this figure are the 
characteristic times and the light-curve scaling laws. 



5 CONCLUSION AND DISCUSSION 

Based on the shock model which has been essentially suc- 
cessful to explain long GRB afterglows, we here derive the 
light curves of short GRB afterglows in the early phase when 
the blast wave is not decelerated by the ambient medium 
considerably. The reverse-shock emission has been ignored at 
the beginning since it is always Newtonian initially for short 
GRBs. We consider the pair-loading effects on the emission. 
The spectrum shows rapid soft-to-hard evolution at the first 
several seconds (t < tf), and then a usual hard- to-soft evolu- 
tion after several tens of seconds (t > tdcc). Simultaneously, 
there are two peaks appearing in the light curves in the op- 
tical to hard X-ray range. The first "pair peak" will appear 
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Figure 1. Example of early afterglows of from short GRBs at 
two fixed frequencies, v = 2 X 10 14 (upper frame) and 5 X 10 17 
(bottom frame) Hz. The parameter values taken to calculate the 
light curves are: E 1 = = 10 52 ergs, n = 0.01 cm -3 , fa = 
10 2 , p = 2 and the others equal to typical values of long GRBs 
(see details in the text). The characteristic times and the scaling 
laws of fluxes with time are marked. The dashed lines show the 
sensitivity of Swift instruments, the X-ray (XRT) and UV optical 
(UVOT) Telescopes. 



at the optical, but it is too faint (mag ~ 21) and too short 
(~ 2 s) to be detected by any current and upcoming instru- 
ment. But the double-peak structure in the light curve is 
expected to be observed at X-ray band: the first peak at t p k 
and the second peak at around tdcc - It took only 20 to 70 s for 
Swift to point its Narrow Fields instruments, consisting of X- 
ray and UV optical Telescopes, to the GRB direction, short 
GRBs will be easily detected before the second peak (Fig. 
1). The recently proposed micro-satellite ECLAIR (Barret 
2003) is even expected capable of detecting the first peak. 

Though the reverse shock becomes mildly relativistic 
finally at t dcc , we have neglected its emission here, which 
is mainly in the soft band, say, the optical. If we consider 
further the effect of pair-loading in the fireball which due 
to 7 — 7 absorption of the prompt burst in the fireball, it 
would be in much softer band such as in the IR. Because 
the same energy may be shared by more produced e pairs, 
the lower-energy leptons would radiate at softer frequency. 
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So the reverse-shock emission will not affect the X-ray light- 
curve, though may affect the optical one. 

The blast wave emission in sub-MeV must be under 
the BATSE detection threshold for short GRBs. Provided 
that the energy is Eh = 10 52 ergs, and e e ,-i = da, 28 = 1 
similar to typical values of long GRBs, we find a constraint 
on initial Lorentz factor and ambient density of r^^n 1 ^ < ~ 
1 (cf.eq.[17]), and that the Lorentz factor of short GRBs 
is not allowed to be large, i.e., r\ < 10 3 . This also limit 
the ambient density to n 5s 0.1 cm~ 3 , which is consistent 
with upper limit on late-time short GRB afterglows. So far, 
the best constraint on short GRB afterglows comes from 
the observation of short/hard GRB 020531, which yields the 
limiting magnitudes in R band: 18.5 at 88 min and 25.2 
at 2.97 d (Klotz ct al. 2002). Under a standard afterglow 
model, these data do not allow for a dense medium, i.e. 
n ^ 0.1 cm -3 (see Fig. 2 in Panaitescu et al. 2001). Low 
densities favor the GRB model related to compact object 
mergers (Eichler et al. 1989; Narayan, Paczyhski & Piran 
1992) in galactic haloes or in the intergalactic medium. 

At tioad(> 2 s) the optical photons in the pulse may be 
up-scattered to MeV by synchrotron self-Compton process. 
But the flux is of orders lower than BATSE detection thresh- 
old, and is unable to change the short-duration property of 
short GRBs. 

Unlike the long GRBs which may overlap the early af- 
terglows and lead to complication, the short GRBs stop 
abruptly. And due to lower ambient density the blast waves 
take longer time to begin decelerating considerably, so their 
early afterglows are easy to be observed. If detected and con- 
firmed, the double-peak structure in early afterglows has an 
important indication for short GRBs — with redshift hav- 
ing been measured, we can determine the most important 
parameter r\ from eqs.(4) and (11) (Beloborodov 2002), and 
then we can further use the value of rj to constraint the other 
parameters, Ek and n, by eq.(l). So an observation of early 
afterglows provides important constraints on the short GRB 
parameters, related both to the relativistic flow and to the 
surroundings. 
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